An entropy— driven noise— induced phase transition 
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We introduce a class of exactly solvable models which exhibit an ordering noise-induced phase 
transition driven by an entropic mechanism. In contrast with previous studies, order does not 
appear in this case as a result of an instability of the disordered phase produced by noise, but of 
a balance between the relaxing deterministic dynamics and the randomizing (entropic) character 
of the fluctuations. A finite-size scaling analysis of the phase transition reveals that it belongs to 
the universality class of the equilibrium Ising model. All these results are analyzed in the light of 
the nonequilibrium probability distribution of the system, which can be derived analytically. The 
relevance of our result as a possible scenario of the so-called Lower Critical Solution Temperature 
(LCST) transitions is analysed. 
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A surprising discovery in the field of stochastic pro- 
cesses is that the interaction between noise and spatial 
degrees of freedom in extended nonlinear systems can 
lead, under appropriate conditions, to noise-induced or- 
der Q. This seemingly counterintuitive phenomenon 
arises when a distributed medium develops some kind 
of spatiotemporal self-organization which requires noise 
to exist. The resulting ordered state can take many 
forms, such as a static spatial pattern traveling ex- 
citable waves [pi, propagating pulses spatiotemporal 
chaotic states B, or a coherent response to an external 
periodic forcing |(|, among others. From a fundamen- 
tal point of view, however, the most interesting exam- 
ple is that of noise-induced ordering phase transitions 
(NIOPTs) §. In this case, transitions between true ex- 
tended phases (in a thermodynamic sense) are produced 
when the noise intensity is used as control parameter. 
These non-equilibrium phase transitions can therefore be 
characterized using standard techniques of critical phe- 
nomena, such as dynamic renormalization group, mean- 
field approximations, as well as finite-size scaling analy- 
ses using, for instance, results obtained from numerical 
simulations of stochastic partial differential equations. 
Up to now, all the arguments presented to account for 
the occurrence of NIOPTs have been dynamical ones. 
This is mainly because it has been impossible so far to 
find a model whose nonequilibrium steady-state proba- 
bility distribution and the associated effective potential 
are known. In particular, noise-induced phase transitions 
have been systematically explained in terms of a short- 
time instability of the local dynamics, which becomes 
"frozen" at larger times by the spatial coupling (t]|| . 

In this Letter we introduce a class of systems exhibit- 
ing NIOPTs for which the steady-state probability dis- 
tribution can be obtained exactly, so that one can define 
the corresponding nonequilibrium free energy or poten- 
tial. As a consequence, the NIOPT can be studied in the 
steady state, with no dynamical reference. It turns out 



that the noise-induced phase transition that appears in 
these model systems is not a consequence of an instability 
of the disordered homogeneous state. Rather, the situa- 
tion is similar to what happens for noise-induced tran- 
sitions in O-dimensional systems, where the disordered 
state is linearly stable, but the effective non-equilibrium 
potential in the steady state is bimodal |) 10 1. A phase 
transition is obtained by coupling adequately many of 
those O-dimensional systems. For a strong enough spa- 
tial coupling, the neighboring variables tend to a com- 
mon value and a macroscopic ordered phase appears as 
a result of ergodicity breaking. This situation contrasts 
vividly with NIOPTs reported so far, where the need of 
a short-time dynamical instability prevents systems that 
undergo noise-induced transitions in 0-d from exhibiting 
NIOPTs in the presence of spatial coupling 0. 

The phase transition that we report in what follows 
shares some features with the so-called Lower Critical 
Solution Temperature (LCST) transition which appears 
in some polymer blends jnj . In these very special sys- 
tems, one can go from a disordered phase (homogeneous 
mixture) to an ordered phase (separated phases) by in- 
creasing the temperature (inverted phase diagram). This 
is an equilibrium phenomenology which goes in the op- 
posite direction of what is expected intuitively from sta- 
tistical physics. Our findings present a possible scenario, 
in terms of a Langevin description, in which these tran- 
sitions can take place. To reach our goal we consider 
two important theoretical ingredients: the presence of a 
field-dependent kinetic coefficient and the existence of a 
fluctuation-dissipation relation accordingly. 

Let us consider the following generic deterministic 
model for a real field <f>(r, t) : 



dt 
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which corresponds to a relaxational flow [Q in a po- 
tential U({0}) with a field-dependent kinetic coeffi- 
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cient r(0). Field-dependent kinetic coefficients can ap- 
pear in coarse-grained, master-equation derivations of 
macroscopic (deterministic or stochastic) field equations 
[^3|,0. Their explicit form depends on the assumptions 
in their derivation, but the general trend is that the co- 
efficient is large in the disordered phase and small in the 
ordered phase Jf4|jl]|. These coefficients have been in- 
troduced in the context of phase separation in binary 
mixtures, either to explore ad hoc their influence, or by 
necessity to take into account the effect of a gravitational 
field |16|. In those cases, the kinetic coefficients are also 
a function of the spatial derivatives of the field, in or- 
der to reflect the conserved nature of the concentration 
field. Here, on the other hand, we consider for simplicity 
a non-conserved situation. 

When stochastic terms are introduced in this descrip- 
tion a fluctuation-dissipation relation is used jf4|]. Fol- 
lowing this approach we define that the Langevin equa- 
tion (in the Stratonovich interpretation) corresponding 
to the deterministic equation ([l]) to be: 

-_- m r,t))j^+mm 1/2 mt)- (2) 



r(0)-i/(i + c ^ 2 ), 



(7) 



dt 



The noise is taken to be Gaussian, with zero mean and 
correlation 



W,mr',t')) = 2e5(f-f')5{t-t'). 



(3) 



where e is the noise intensity. Under these conditions, 
the stationary solution of the probability distribution for 
the field, P st ({</>}) is of the Boltzmann's type: 



PM}) ~ e - y -/ £ 

in terms of an effective potential 



V cS (W) = V ({<!>}) + j J dflnT^r)) , 



(4) 



(5) 



where £o is a renormalized parameter, proportional to e, 
which includes an ultraviolet cutoff |17| ], The previous 
expression defines the exact nonequilibrium free energy 
of model (||). We will now show that this class of systems 
can exhibit a noise-induced phase transition. 

So far, the discussion has been made in very general 
terms. In order to proceed further, we assume that the 
deterministic free-energy potential has the usual Landau 
form 



dr 



V {<l>{f)) 



K 



V</>(r)| : 



(6) 



and we choose a monostable local potential V~o{4>) = § <j> 2 , 
with a > 0. We also consider a non-conserved order pa- 
rameter system in which the kinetic coefficient T is only a 
function of the field and not of its spatial derivatives. Ac- 
cording to the discussion above, we adopt the functional 
form 



which means that fluctuations are larger in dilute regions 
and smaller in dense ones. In this particular case, and 
according to Eqs. (^) and (||), the local effective potential 
can be seen to become bistable for £o > a/c. One can 
then expect that for a strong enough coupling K, a true 
phase transition towards an ordered state controlled by 
the noise intensity e exists. In order to check this expec- 
tation, we perform a Weiss mean-field analysis of the cor- 
responding local effective potential fll8| . The analysis is 
more easily done by considering a spatial discretization in 
a regular d-dimensional lattice with spacing Ax = 1, i.e. 
4>(fi) — > fa, with i the cell index. The gradient term is 
then approximated by the sum over nearest neighbors on 

the lattice in the standard way, |V</>| 2 — > J2j=i (<Aj — 4>i) 2 - 
The mean-field approximation is tantamount to replacing 
the exact value of the neighbors by a mean-field common 
value M. In this way, the local effective potential be- 
comes 

V eS (<j>, M) = V (<P) + y lnr(^) +dK{4>~ Mf , (8) 

and depends on the yet unknown mean-field value M. 
This can be computed by using the consistency relation: 



M = 



J — oo t 



-Kff(0,M)/e 



I: 



-VM<t>,M)/e c 



(9) 



which, in general, has to be solved numerically. The so- 
lution of this equation is shown as a solid line in Fig. |l| 
for the two-dimensional case. The result predicts a phase 
transition from disorder (M = 0) to order (M ^ 0) as 
the noise intensity increases. This theoretical result has 
been confirmed by numerical simulations of a discretized 
version of model (|^) in a two-dimensional square lattice. 
Following standard procedures of equilibrium statistical 
mechanics, we have characterized the phase transition by 
means of two quantities: the intensive order parameter 
M = (m), where m — \ 'Yl li 4>i\/ L 2 is the absolute value 
of the spatial average of the field, and the generalized 
susceptibility 



X = L' 



(to 2 ) — (m) 2 



(10) 



which measures the relative fluctuations of the field in 
its steady state (e plays the role of KbT in equilibrium 
phase transitions and L is the linear system size). The 
numerical results are plotted in Figs. |f| and ||, which show 
the behavior of the order parameter M and the suscep- 
tibility x, respectively, as the noise intensity e increases, 
for four different system sizes. The figures clearly denote 
the existence of a noise-induced ordering phase transi- 
tion with the usual features of a second-order equilib- 
rium phase transition, in agreement with the theoretical 
analysis. Monte Carlo simulations of the system making 
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use of the nonequilibrium effective free energy give com- 
pletely compatible results. The position of the critical 
point has been computed both by extrapolating to infi- 
nite size the location of the maxima of the susceptibility 
and by means of a finite-size scaling analysis (see below) . 
The two approaches lead to similar results. 
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FIG. 1. Intensive order parameter M as function of the 
noise intensity for o = 1, c = 0.5, K = 1 and for four differ- 
ent system sizes. The continuous line is the Weiss mean-field 
result. 
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FIG. 2. Generalized susceptibility \ versus noise intensity 
for the case of Fig. [I]. 

As mentioned above, noise-induced phase transitions 
have been related so far to an instability of the homoge- 
neous phase at short times. This had been clearly mani- 
fested by a linear stability analysis of the equation for the 
first moment, and was understood as a renormalization 
of the linear forces due to the noise term. Let us now 
perform in our model a similar linear stability analysis of 
the homogeneous disordered state = jl9],|o) and we 
will be able to show that this is not the mechanism that 
induces the phase transition in our model. The dynami- 
cal equation for the first statistical moment, in the linear 
approximation, reads 



at 



= -(a + e c)(0)+irV 2 (</») 



(11) 



The contribution of the noise is in the same direction 
as the one from the local potential, and thus the homo- 
geneous state 4> — is linearly stable. Hence the phys- 
ical mechanism responsible for the appearance of a new 
phase must be different from the one of previously known 
cases of noise induced phase transitions. A simple inter- 
pretation can be found in terms of a balance between 
the role of the deterministic monostable local potential, 
which tends to take the system towards the disordered 
(4> = 0) phase, and the stochastic motion, which is due 
to the fact that fluctuations are more intense in the dis- 
ordered phase, and thus push the system away from it. 
The presence of a spatial coupling K helps to break the 
symmetry of the homogeneous state. 

A natural question arises concerning the universal- 
ity class of the noise-induced phase transition described 
above. The availability of an exact nonequilibrium po- 
tential for this model, given by Eq. (|J), allows us in a 
very simple way to analyze its universal behavior at the 
critical point. If an expansion of In T (</)), as chosen in Eq. 
(0), is performed on the nonequilibrium effective po- 
tential can be seen to have a Ginzburg-Landau form with 
infinite higher-order (irrelevant) terms. Since the univer- 
sality class of the Ginzburg-Landau model is that of the 
Ising model, we predict that this transition will belong to 
the d-dimensional Ising universality class. This is indeed 
confirmed by a finite-size scaling analysis of the numer- 
ical results. We make use of the following scaling laws, 
generalized from equilibrium statistical mechanics JD: 



M = L- f} ' v M{el}/ V ), 



(12) 
(13) 



where Novikov's theorem has been applied to the term 



where M and x are the corresponding scaling functions, 
and £ is the reduced control parameter, e = (1 — e/e c ). 
The values of the critical exponents /3, 7 and v define 
the universality class. We have chosen the correspond- 
ing values of the 2-d Ising universality class, and plotted 
expressions ( |l2|) and (|l3|) in Figs. || and f|. All data are 
seen to scale fairly well into the same curve in the region 
close to the critical point. As expected, those data corre- 
sponding to the smallest system size (L = 16) show the 
major discrepancies. Hence, we can conclude that this 
new kind of noise-induced phase transitions belongs to 
the 2d-Ising universality class. 

In conclusion, we have introduced a dynamical system 
that exhibits a nonequilibrium phase transition from dis- 
order to order as the intensity of the noise is increased. As 
a consequence, the system is disordered below a certain 
amount of noise and can be considered a possible sce- 
nario, within the context of stochastic differential equa- 
tions, for the appearance of a lower critical tempera- 
ture in physical systems. At variance with other noise- 
induced nonequilibrium phase transitions reported be- 
fore, the physical mechanism inducing the transition is 
not a noise-induced short time instability of the disor- 
dered state. Rather, we are able to show that the effec- 
tive local potential, which can be obtained analytically, 
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exhibits bistability as the intensity of the noise increases, 
whereas the deterministic one is monostable. Hence, the 
mechanism is seen to be given by a balance between 
deterministic forces and stochastic motion. Our model 
shows that it is possible to obtain a noise-induced phase 
transition by coupling adequately O-dimensional systems 
undergoing a noise-induced transition. We have argued 
that this type of coupling is consistent with a mesoscopic 
description in terms of fields, kinetic coefficients and the 
fluctuation-dissipation relation. Our findings have been 
sustained by numerical simulations in a two-dimensional 
lattice model. Finally, by studying the exact stationary 
probability distribution, we are able to show that the 
model belongs to the same universality class as the Ising 
model. 
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FIG. 3. Finite-size scaling function of the order parameter 
M with the 2d-Ising critical exponents. 
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FIG. 4. Finite-size scaling function of the susceptibility \ 
with the 2d-Ising critical exponents. 
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